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Abstract 

The recent detection of a hard X-ray component in the supernova remnant Cassiopeia A is interpreted as 
synchrotron emission from electrons accelerated to energies up to 40 TeV (Allen et al., 1997). It is therefore 
tempting to consider TeV gamma-ray emission from this object through : i) bremsstrahlung and inverse Comp- 
ton radiation from electrons and/or ii) vr'^ production from an associated high energy cosmic ray component 
hitting surrounding material. Cas A was observed by the CAT imaging Cherenkov telescope during the ob- 
serving season Aug-Nov 1998. An upper limit to the integral flux above 400 GeV of 0.74 x 10~^^ 7cm^^s^^ 

0^ ■ is derived. This result is used to constrain shock-acceleration models for production of VHE gamma-rays in 

^ ■ SNRs. 

j^ '- 1 Introduction: 

^ . The supernova remnant Cassiopeia A (Cas A), aged ~300 yr at a distance of -^2.8 kpc, is the brighest radio 

source in the sky. The radio emission region is concentrated in a shell of radius ~130 arcsec and thickness 
~30 arcsec. X-ray data from the RXTE satellite indicate the presence of electrons accelerated to energies up to 
T-H ■ ~40TeV which emit synchrotron photons above 10 keV (Allen et al., 1997). These electrons could generate 7- 

^ \ rays in the GeV-TeV band via inverse Compton or bremsstrahlung mechanisms. The 7-ray emissivity critically 

^ I depends on the value of the magnetic field in the acceleration region. In addition, a nuclear component is likely 

O ' to be accelerated along with electrons giving rise to additional 7-rays produced in nuclear collisions via vr'^ 

p^ , decay. Upper limits to the Cas A 7-ray flux above 100 MeV were obtained by the satellite experiments SAS-II 

Q>^ \ and COS-B. Analyzing these data, Cowsik and Sarkar (1980) were able to set a lower limit of 8 x 10~^G to the 

<3\ • magnetic field in Cas A. Further observations by the EGRET experiment on the CGRO satellite led Esposito 

^ . et al. (1996) to derive a more constraining lower limit of 2x 10~^G. 

Qh| The situation in Cas A is reminiscent of that in SN1006 where Koyama et al. (1995) inferred the presence 

O ' of up to 100 TeV electrons from the observed X-ray synchrotron emission. Subsequently, TeV 7-ray emission 

-M ■ was observed by the CANGAROO atmospheric Cherenkov telescope (Tanimori et al., 1998). The Whipple 

j^ \ collaboration, analyzing a 2.7 hr exposure on Cas A, reported an upper limit to the 7-ray flux above 300 GeV 

l~^ ' of 2.6xl0^^^cm^^s^^ (Lessard, 1996). More recently, the same group derived an improved upper limit of 

1.8xlO~^^cm~^s^^ above 400 GeV from a 6.9 hr exposure (Lessard, 1999). In this paper, we present an 
analysis of a 24.4 hr exposure on Cas A obtained with the imaging atmospheric Cherenkov telescope CAT. 
C^ I The CAT telescope and the data analysis technique are described in section 2. The database is presented in 

section 3. The results of the analysis are summarized in section 4. The implications for current models of 
7-ray production in SNR's are discussed in section 5. 

2 The CAT telescope and data analysis: 

The CAT imaging atmospheric Cherenkov telescope is described in Barrau et al. (1998). It consists of a 
18m^ Davies-Cotton mirror fitted with a 600 PMT fine-grained camera (pixel size: 0.12° diameter) at the focal 
plane. The energy threshold is 250 GeV at zenith. The telescope is located in the French Pyrenees at an altitude 
of 1650m and geographical coordinates 2°E, 42°N. The data analysis technique used in the present paper is 
detailed in Le Bohec et al. (1998). The basis for the analysis of individual shower images is to compare the 
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observed image to a template oi calculated images at diiierent energies, zemth angles and impact parameters 
to the telescope. Gamma/hadron discrimination is performed with a cut on a single x^-like parameter. The 
analysis also yields the energy, impact parameter and orientation (the equivalent of the a-parameter in the 
standard Hillas method) for each event. Note that, in the present study, Cas A is considered as a point-like 
source, just as the Crab nebula is. Extensive Monte-Carlo simulations of the instrument have been performed 
and proved highly reliable when checked against real 7-rays from the Crab or Mkn501 (e.g. Goret, 1997, 
Mohanty, 1999). 

3 The Cas A database: 

The observations relevant for the present analysis were performed during the period August-November 
1998. ON-source runs lasted for 30 minutes followed or preceded (but not systematically) by OFF-source 
runs shifted by it 35 minutes in RA. After selection on meteorological conditions, zenith angle (<25°) and 
correct experiment operation, a total exposure of 24.4 hr ON-source and 13.6 hr OFF-source were retained for 
further analysis. A database on the Crab nebula, with the same selection criteria and observing period, was 
also compiled to allow for cross-checking of the analysis. The trigger condition required that 4 PMTs above a 
threshold of 2.5 photoelectrons fire simultaneously within a sector (see Barrau et al., (1998), for details). The 
average raw trigger rate for these two data samples was --^17 Hz. Both the Cas A and Crab data were processed 
the same way as is decribed in the next section. 
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Figure 1 : The distributions of the orientation angle a for the Crab (a) and Cas A (b) datasets respectively 
using the x^ method. The ON and (renormalized) OFF-source data are represented by the solid histogram 
and stars respectively. The insets shows the ON-OFF distributions. The data are binned in 2.5° intervals. 



4 Results: 

The Crab data were used to optimize the cuts for the best sensitivity in the search for a signal from Cas A. 
This analysis yielded the following cuts for gamma/hadron (shape) discrimination: i) x^ probability >0.40 and 
ii) a total number of photoelectrons in the image of Qtot>50 p.e.. In addition, the estimated impact parameter 
was restricted to <160 m. The a-plots for the Crab and Cas A are displayed in figure 1. The OFF data were 



normalized using the tracJiing ratio, i.e., the ratio noN/noFF oi the number oi events surviving the shape cuts 
in the region 20° < a <120°. The orientation cut was set to amax=^°, which retains --^ 80% of the 7-rays 
selected after the shape cuts according to Monte-Carlo simulations. The results of the analysis are summarized 
in Table 1 . 





Crab 


CasA 


Ton (hrs) 


13.3 


24.4 


Toff (hrs) 


9.2 


13.6 


NoN (events) 


2006 


2360 


NoFF (events) 


706 


1259 


Tracking ratio 


1.43 


1.85 


Excess (events) 


998 


33.5 


Standard Deviation (events) 


60.7 


85.9 


Significance (n^) 


16.45 


0.39 



Table 1. Results of the analysis for the Crab and Cas A. 

^From these data, a 3cr upper limit to the Cas A flux can be set at <0.14 times the Crab flux. Monte-Carlo 
calculations indicate that the energy threshold subsequent to the above cuts is 400 GeV, with ~89% of the 
selected events in the sample having an effective energy above 400 GeV. The effective area, assuming a Crab- 
like energy spectrum for Cas A, is estimated to be ~3.5xlO®cm^. With the assumption of a null detected 
signal, this translates into a 3o- upper limit to the Cas A 7-ray flux of: 

$(>400 GeV) < 0.74xlO~^Scm~2s-i 

The implications of this result concerning the current models for TeV-gamma production in SNRs are exam- 
ined in the next section. 

5 Discussion: 

As was already noted by Cowsik and Sarkar (1980), considering 7-ray production via bremsstralhung 
and inverse Compton, upper limits on the 7-ray flux from Cas A must constrain primarily the strength of the 
magnetic field in the source. The problem was recently reconsidered by Baring et al. (1999) using a non-linear 
shock-acceleration model. The model aims at predicting the overall continuum photon emission spectrum in 
SNRs from radio to 7-rays. The model was tailored to the specific case of Cas A by Ellison et al (1999), 
taking into account relevant observations over the full photon energy range. The results are displayed in figure 
2, which shows the best fit of the model to the data from radio through X-rays up to 7-rays. The contributions 
of the different components (synchrotron, bremsstrahlung, IC and p-p) in each energy band are indicated. The 
conclusion is that in order to account for the upper limits reported in both the GeV and TeV ranges, the strength 
of the magnetic field must be greater than ~10^^G. As discussed by Cowsik and Sarkar (1980), such a high 
value in excess of the equipartition field of '-^4x lO^^G calls for a magnetohydrodynamic field amplification 
which remains to be fully investigated (see Keohane, (1998), for a recent discussion). The present result may 
also be considered in the frame of the diffusive shock-acceleration model of Drury, Aharonian and Volk (1994) 
as a test of cosmic-ray acceleration in SNRs. 

6 Conclusions: 

The improved upper limit to the TeV 7-ray flux from the SNR Cas A presented in this paper sets a strin- 
gent lower limit to the magnetic field in the acceleration region. This result should help understanding the 
mechanisms at work in this young object. 
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Figure 2: Predictions of the nonlinear diffusive shock-acceleration model from radio to TeV ^-rays as com- 
pared to observations (see Ellison et al. (1 999) for details and references therein). The present result is shown 
together with the Whipple upper limit ofLessard (1999). 
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